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Beschreibung 

Verfahren zum Herstellen eines Volumendatensat zes 

5 Die Erfindung betrifft ein Verfahren zum Herstellen eines Vo- 
lumendatensat zes . 

Insbesondere mit modernen bildgebenden medizintechnischen Ge- 
raten aufgenoramene Bilder weisen eine relativ hohe Auflosung 

10 auf, so dass mit ihnen verstarkt 3D-Aufnahmen (Volumendaten- 
satze) erstellt werden. Bildgebende medizintechnische Gerate 
^ sind z.B. Ultraschall-, Computertomographie- , Magnetresonanz- 
oder Rontgengerate oder PET-Scanner. Ferner konnen ofter Com- 
putertomographie- (CT) oder Rontgengerate eingesetzt werden, 

15 da sich eine Strahlenbelastung, die ein Lebewesen wahrend ei- 
ner Untersuchung mit einem dieser Geraten ausgesetzt ist, 
verringert hat. Volumendatensat ze weisen jedoch eine groftere 
Datenmenge auf als Bilddatensat ze von herkommlichen zweidi- 
mensionalen Bildern, weshalb eine Auswertung von Volumenda- 

20 tensatzen relativ zeitauf wandig ist. Die eigentliche Aufnahme 
der Volumendatensat ze dauert zur Zeit in etwa eine halbe Mi- 
nute, wobei man fur das Durchforsten und Aufbereiten des Vo- 
lumendatensat zes oft eine halbe Stunde oder mehr benotigt. 
Daher sind automatische Erkennungs- und Auf bereitungsverf ah- 
^£J5 ren notwendig und willkommen. 

Bis etwa zum Jahr 2000 war es in der Computertomographie (CT) 
fast nur iiblich, eine Diagnose anhand axialer Schichtstapel 
(Schnittbilder) zu treffen oder sich zumindest fur einen Be- 

30 fund vorwiegend an den Schnittbildern zu orientieren. Seit 
etwa 1995 sind dank der Rechenleistung von Computern 3D- 
Darstellungen auf Bef undungskonsolen verbreitet; sie hatten 
aber zuerst eher wissenschaf tliche oder erganzende Bedeutung. 
Um dem Arzt eine Diagnose zu erleichtern, sind ferner im We- 

35 sentlichen vier Grundverf ahren der 3D-Visualisierung entwi- 
ckelt worden: 
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1. Multiplanare Ref ormatierung (MPR) : Dies ist nichts anderes 
als eine Neuzusammenstellung des Volumendatensat zes in ande- 
rer Orientierung als z.B. den urspriinglichen horizontalen 
Schichten. Es wird insbesondere zwischen der orthogonalen MPR 

.5 (3 MPRs , jeweils senkrecht zu einer Koordinatenachse) , der 

freien MPR (schrage Schichten; abgeleitet = interpoliert ) und 
der Curved MPR ( Schichterstellung parallel zu einem beliebi- 
gen Pfad durch das Abbild des Korpers des Lebewesens und z.B. 
senkrecht zu der MPR, in welcher der Pfad gezeichnet wurde) 

10 unterschieden . 



2. Shaded Surface Display (SSD) : Segmentierung des Volumenda- 
tensatzes und Darstellung der Oberflache der herausgeschnit- 
tenen Objekte, meist stark gepragt durch Orientierung an den 

15 Grauwerten der Abbildung (z.B. CT-Werten) und manuelles 
Hilf s-Editing. 

3. Maximal Intensity Projection (MIP) : Darstellung der hochs- 
ten Intensitat entlang jedes Sehstrahls. Bei der so genannten 

20 Thin MIP wird nur ein Teilvolumen dargestellt. 

4. Volume Rendering (VR) : Darunter wird eine Modellierung der 
Abschwachung des Sehstrahls, der vergleichbar einem Rontgen- 
strahl in das Objekt eindringt, verstanden. Dadurch wird die 

^kL5 gesamte Tiefe des abgebildeten Korpers (teilweise durchschei- 
nend) erfasst; es gehen jedoch Einzelheiten von kleinen und 
vor Allem dlinnschichtig dargestellten Objekten verloren. Die 
Darstellung wird manuell durch Einstellung so genannter 
Transf erf unktionen ( Farb-Lookup-Tabellen) gepragt . 

30 

Obwohl das Volume Rendering (VR) relativ gute Resultate bie- 
tet, konnen mit dem Volume Rendering bearbeitete Bilder eine 
begrenzte Transparenz aufweisen. Durch eine Veranderung der 
Transf erf unktion kann namlich die Transparenz nur mit Ein- 
35 schrankungen iiber den gesamten Bereich gewahrleistet werden. 
Es kann vorkommen, dass die alleinige transparente Darstel- 
lung bestimmter erkrankter Strukturen, wie z.B. ein erkrank- 
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tes Organ, nur begrenzt durch eine Einstellung der Transfer- 
funktion moglich ist. Dann konnen Editierf unktionen im Umfeld 
des Volume Rendering herangezogen werden, um zumindest mit 
manueller Unterstutzung die gewunschten abgebildeten Struktu- 
5 ren herauszusegmentieren und dann separat dar zustellen . Oft 
mlissen die gewunschten abgebildeten Strukturen komplett manu- 
ell segmentiert werden. Sind die Strukturen mit einem Volu- 
mendatensatz, der in Form mehrerer auf einanderf olgender com- 
putertomographischer Schnittbilder vorliegt, abgebildet, so 
10 werden die Konturen der abgebildeten Strukturen Schicht fur 
Schicht manuell umfahren, so dass nur die in den geschlosse- 
nen Konturen enthaltenen Voxel dargestellt werden. 

Die Aufgabe der Erfindung ist es daher, ein Verfahren zum 
15 Herstellen eines Volumendatensat zes aus einem ein abgebilde- 
tes Objekt umfassenden Volumendatensat z anzugeben, bei dem 
interessierende Strukturen des abgebildeten Objekts verbes- 
sert dargestellt werden. 

20 Die Aufgabe der Erfindung wird gelost durch ein Verfahren zum 
Herstellen eines Volumendatensat zes , aufweisend folgende Ver- 
f ahrensschritte : 

- Segmentieren der abgebildeten Oberflache eines in einem 
ersten Volumendatensat z abgebildeten ersten Objektes, 

- Transf ormieren des ersten Volumendatensat zes in einen zwei- 
ten Volumendatensat z derart, dass die segmentierte abgebil- 
dete Oberflache in eine Ebene transf ormiert wird, 

30 

- Herstellen eines dritten Volumendatensat zes, indem der 
zweite Volumendatensat z derart gefiltert wird, dass nicht- 
interessierende, im zweiten Volumendatensat z abgebildete 
Strukturen des ersten Objektes aufgrund nichtinteressieren- 

35 den Strukturen zugeordneter Merkmale und aufgrund der von 

der Oberflache zu erwartenden Entfernungen der nichtinte- 
ressierenden Strukturen herausgef iltert werden und interes- 



200214682 



4 

sierende, im zweiten Volumendatensatz abgebildete Struktu- 
ren des abgebildeten ersten Objektes aufgrund interessie- 
renden Strukturen zugeordneter Merkmaler und aufgrund den 
von der Oberflache zu erwartenden Entfernungen der interes- 
5 sierenden Strukturen erhalten bleiben. 

Im ersten Volumendatensatz ist zumindest ein Teil des ersten 
Objektes, das beispielsweise ein Lebewesen ist, abgebildet. 
Sollen nun mit einem dem ersten Volumendatensatz zugeordneten 

10 Bild weiter im Inneren des ersten Objektes liegende Struktu- 
ren untersucht werden, so werfen die Abbilder weiter an der 
^ Oberflache angeordneter Strukturen des ersten Objektes ihre 
Schatten auf die weiter innen abgebildeten Strukturen oder 
verschleiern deren Darstellung. Mittels des erf indungsgemalien 

15 Verfahrens sollen moglichst die weiter auften angeordneten ab- 
gebildeten Strukturen herausgef iltert werden, ohne abgebilde- 
te Strukturen, die untersucht werden sollen und weiter im In- 
neren des ersten Objektes angeordnet sind, zu entfernen. 

20 Zunachst wird erf indungsgemafi insbesondere automatisch die im 
ersten Volumendatensatz abgebildete Oberflache des abgebilde- 
ten ersten Objektes bestimmt (heraussegmentiert ) . Im Falle 
eines Lebewesens als erstes Objekt wird diese in der Regel 
gekrummte Oberflache nun in die Ebene transf ormiert , als ob 
v ^£5 man das abgebildete erste Objekt abrollen wurde. Man denke 
' hier zum Vergleich an die Projektion der Erdoberf lache auf 

Landkarten. Insbesondere wenn es sich bei dem abgebildeten 
ersten Objekt urn den abgebildeten Torso des Lebewesens han- 
delt, der quasi Saulenform, mit naherungsweise elliptischer 

30 Grundflache aufweist, lasst sich die abgebildete Oberflache, 
d.h. die abgebildete Korperoberf lache des Lebewesens in eine 
ebene Flache abrollen. Somit entsteht der zweite Volumenda- 
tensatz. 

35 Anschlieftend werden die herauszuf ilternden Strukturen mit ei- 
nem geeigneten Filter oder einem geeigneten Set von Filtern 
aus dem zweiten Volumendatensatz herausgef iltert . So wird fur 
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jede Art oder Klasse von herauszuf ilternden Strukturen, die 
beispielsweise Haut, Fett, Rippen, Knochen oder Muskeln sind, 
ein Filter ermittelt. Die jeweiligen Filter werden nach den 
herauszuf ilternden Strukturen zugeordneten Merkmalen entwi- 
5 ckelt. Dabei wird gemali einer Ausf iihrungsf orm der Erfindung 
den nichtinteressierenden Strukturen zugeordnete dichteorien- 
tierte, texturorientierte, kantenbetonte und/oder morphologi- 
sche Filterung verwendet. Die einzelnen Filterantworten kon- 
nen letztlich in einer Mer kmalsmatrix geeignet zusammenge- 
10 rechnet werden. 

^ Des Weiteren sollen Strukturen weiter im Inneren des ersten 
Objektes untersucht werden. Daher wird fur die Filterung zum 
Herausf iltern der nichtinteressierenden Strukturen zusatzlich 
15 die erwarteten Entfernungen der herauszuf ilternden Strukturen 
von der Oberf lache des ersten Objektes berucksichtigt . 

Handelt es sich bei dem ersten Objekt um ein Lebewesen und 
umfasst der Teil des abgebildeten Korpers z.B. den Bauchraum 
20 eines Menschen und die interessierenden Strukturen sind die 
Wirbelsaule und die inneren Organe des Menschen, so umfassen 
nichtinteressierende Strukturen Haut, Rippen und Fettgewebe. 
Fur das Herausf iltern der abgebildeten Haut, der abgebildeten 
Rippen und des abgebildeten Fettgewebes wird also jeweils ein 
{•25 Filter bestimmt, der das Abbild der entsprechenden Struktur 
erkennt . 

Bei dem Filter, der die abgebildeten Rippen herausf iltern 
soli, handelt es sich beispielsweise um einen Filter, der ab- 

30 gebildete Knochen erkennen soli. Damit jedoch die abgebildete 
Wirbelsaule, die ein Teil der interessierenden Strukturen 
ist, nicht beeintracht igt wird, soli der abgebildeten Knochen 
zugeordnete Filter nur fur einen Bereich des abgebildeten 
Korpers wirken, der einer bestimmten Entfernung von der Kor- 

35 peroberf lache in das Korperinnere entspricht. Diese Entfer- 
nung entspricht der zu erwartenden Entfernung der Rippen von 
der Korperoberf lache in das Innere des Menschen. Somit ist 
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sicher gestellt, dass kein Teil der abgebildeten Wirbelsaule 
irrtiimlich herausgef iltert wird. 

Der der Haut zugeordnete Filter soli dagegen nur fur einen 
Bereich des abgebildeten Korpers wirken, der der Dicke der 
Haut des Menschen entspricht. Der dem Fettgewebe zugeordnete 
Filter soil ebenfalls nur fur einen Bereich des abgebildeten 
Korpers wirken, der einer vorbestimmten Entfernung von der 
Korperoberf lache in das Korperinnere entspricht. Diese Ent- 
fernung wird derart gewahlt, dass kein Bereich des abgebilde- 
ten Korpers beeintrachtigt wird, in dem interessierende 
Strukturen abgebildet sind. 

Fur die einzelnen Filter wird also u.a. die Topologie des 
Korperinneren des Lebewesens und der zu erwartende Bereich 
der interessierenden Strukturen berucksichtigt . Bei den inte- 
ressierenden Strukturen handelt es sich urn weiter im Korper- 
inneren des Lebewesens angeordnete Strukturen, wie z.B. inne- 
re Organe. Daher sollen oberf lachennahe Strukturen herausge- 
f iltert werden. Fur die Beriicksichtigung der von der Korper- 
oberflache zu erwartenden Entfernungen der nichtinteressie- 
renden Strukturen im Korper des Lebwesens konnen daher bei- 
spielsweise Ortswahrscheinlichkeiten (hoheres Gewicht nahe 
der Oberflache) der nichtinteressierenden Strukturen verwen- 
det werden. 

Die den nichtinteressierenden Strukturen zugeordneten Merkma- 
le werden beispielsweise durch eine Kantendef init ionen in 
Richtung von der abgebildeten Korperoberf lache in das Innere 
des abgebildeten Korpers bestimmt. So ergibt sich z.B. ein 
relativ stark ausgepragter Intensitatsabf all zur abgebildeten 
Lunge oder zu oberf lachennahen Gasblasen im Bauchraum. 

Das erste Objekt kann aber insbesondere auch ein technisches 
Objekt sein, bei dessen Abbild z.B. eine abgebildete Ummante- 
lung oder abgebildete Isolierung des technischen Objektes als 
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nichtinteressierende Strukturen herausgef iltert werden sol- 
len. 

Wenn nach einer Ausf uhrungsf orm der Erfindung der erste Volu- 
5 mendatensatz wenigstens ein abgebildetes zweites Objekt, das 
aulierhalb des ers.ten Objekts angeordnet ist, umfasst, kann 
das abgebildete zweite Objekt aus dem zweiten Volumendaten- 
satz mit den nichtinteressierenden abgebildeten Strukturen 
herausgef iltert werden. Das zweite Objekt ist z.B. ein Tisch, 
10 auf dem das Lebewesen wahrend der Aufnahme des ersten Volu- 
mendatensatzes liegt, Kleidung des Lebewesens oder an dem 
^ Lebwesens angeordnete Instrumente. 

Die abgebildete Oberflache kann, wie es nach einer besonders 
15 bevorzugten Variante der Erfindung vorgesehen ist, vorteil- 
haft segmentiert werden, wenn der erste Volumendatensat z in 
Form mehrerer auf einanderf olgender computertomographischer 
Schnittbilder vorliegt oder als Schichtstapel betrachtet 
wird, die Bilddaten jedes Schnittbildes mit kartesischen Ko- 
20 ordinaten beschrieben sind und fur die Segmentierung der ab- 
gebildeten Korperoberf lache folgende Verf ahrensschritte aus- 
gefiihrt werden: 

- Durchfuhren einer Koordinatentransf ormation fur jedes 
-p|£5 Schnittbild nach Polarkoordinaten bezuglich einer Geraden, 

die durch das abgebildete erste Objekt verlauft und wenigs- 
tens im Wesentlichen rechtwinklig zu den einzelnen Schnitt- 
bildern ausgerichtet ist, 

30 - Ermitteln der Konturen, die in jedem transf ormierten 

Schnittbild abgebildet und der abgebildeten Oberflache zu- 
geordnet sind, 

- Rucktransf ormieren der Bildpunkte der ermittelten Konturen 
35 in das dem ersten Volumendatensat z zugeordnete Koordinaten- 

system und 
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- Re-Extrahieren von Bildpunkten entlang der Konturen fiir die 
Darstellung der in die Ebene transf ormierten Oberflache des 
abgebildeten ersten Objektes. 

GemafJ einer Ausf uhrungsf orm der Erfindung wird aufierdem ein 
vierter Volumendatensat z hergestellt, indem die Bildpunkte 
des dritte Volumendatensat zes in das dem ersten Volumendaten- 
sat z zugeordnete Koordinatensystem rucktransf ormiert werden. 
Der vierte Volumendatensat z umfasst somit die im ersten Volu- 
mendatensat z abgebildeten interessierenden Strukturen . 

Der vierte Volumendatensat z kann nach einer besonders vor- 
teilhaften Variante der Erfindung verwendet werden, urn ein 
dem vierten Volumendatensat z zugeordnetes Bild mittels Volume 
Rendering (VR) dar zustellen . Somit haben nicht herausgefil- 
terte nichtinteressierende Strukturen kaum einen negativen 
Einfluss, weil fiir die eigentlich dar zustellenden Strukturen 
die Transf erf unktion die wesentliche Aussteuerung ist und der 
Charakter der Blending-Technik sehr unempf indlich ist, wenn 
sich "im Vorlauf" jedes Sehstrahls beispielsweise einige we- 
nige Millimeter nichtinteressierende Strukturen befinden. Die 
im vierten Volumendatensat z noch vorhandenen nichtinteressie- 
rende Strukturen sind mit grofter Wahrscheinlichkeit kaum 
sichtbar, da relativ kontrastreiche nichtinteressierende 
Strukturen in der Regel herausgef iltert werden. 

Ein Ausf uhrungsbeispiel ist exemplarisch in den beigefugten 
schematischen Zeichnungen dargestellt. Es zeigen: 

Figur 1 einen Computertomographen, 

Figur 2 einen Volumendatensat z des Bauchraumes eines Patien- 
ten in Form eines aus mehreren Schnittbildern beste- 
henden Volumendatensat zes , 

Figur 3 ein Schnittbild des in der Figur 2 gezeigten Volumen- 
datensatzes, 
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Figur 4 nach Polar koordinaten transf ormierte Bildinf ormatio- 
nen des in der Figur 3 dargestellten Schnittbildes , 

5 Figur 5 einen weiteren Volumendatensatz, bei dem die abgebil- 
dete Korperoberf lache des in der Figur 2 dargestell- 
ten Volumendatensat zes in eine Ebene transf ormiert 
ist, 

10 Figur 6 einen weiteren Volumendatensatz, der moglichst nur 

abgebildete interessierende Strukturen des in der Fi- 
^ gur 5 dargestellten Volumendatensatz umfasst, 

Figur 7 eine mittels Volume Rendering bearbeitete Darstellung 
15 des in der Figur 2 dargestellten Volumendatensat zes , 

und 

Figur 8 eine mittels Volume Rendering bearbeitete Darstellung 
des in der Figur 2 dargestellten Volumendatensat zes , 
20 bei dem abgebildete nichtinteressierende Strukturen 

des in der Figur 2 dargestellten Volumendatensat zes 
herausgef iltert sind. 

Die Figur 1 zeigt schematisch einen Computertomographen mit 
^P^S einer Rontgenstrahlenquelle 1, von dem ein pyramidenf ormiges 
Rontgenstrahlenbundel 2, dessen Randstrahlen in der Figur 1 
strichpunktiert dargestellt sind, ausgeht, das ein Untersu- 
chungsobjekt , beispielsweise einen Patienten 3, durchsetzt 
und auf einen Strahlungsdetektor 4 trifft. Die Rontgenstrah- 

30 lenquelle 1 und der Strahlungsdetektor 4 sind im Falle des 
vorliegenden Ausf uhrungsbeispieles an einer ringformigen 
Gantry 5 einander gegeniiberliegend angeordnet. Die Gantry 5 
ist beziiglich einer Systemachse 6, welche durch den Mittel- 
punkt der ringformigen Gantry 5 verlauft, an einer in der Fi- 

35 gur 1 nicht dargestellten Halterungsvorrichtung drehbar ge- 
lagert (vgl. Pfeil a). 
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Der Patient 3 liegt im Falle des vorliegenden Ausf iihrungsbei- 
spieles auf einem fur Rontgenstrahlung transparenten Tisch 7, 
welcher mittels einer in der Figur 1 ebenfalls nicht darge- 
stellten Tragevorrichtung langs der Systemachse 6 verschieb- 
5 bar gelagert ist (vgl. Pfeil b) . 

Die Rontgenstrahlenquelle 1 und der Strahlungsdetektor 4 bil- 
den somit ein Messsystem, das beziiglich der Systemachse 6 
drehbar und entlang der Systemachse 6 relativ zum Patienten 3 

10 verschiebbar ist, so dass der Patient 3 unter verschiedenen 

Pro j ektionswinkeln und verschiedenen Positionen beziiglich der 

^ Systemachse 6 durchstrahlt werden kann. Aus den dabei auftre- 
tenden Ausgangssignalen des Strahlungsdetektors 4 bildet ein 
Datenerf assungssystem 9 Messwerte, die einem Rechner 11 zuge- 

15 fuhrt werden, der mittels dem Fachmann bekannten Verfahren 
ein Bild des Patienten 3 berechnet, das wiederum auf einem 
mit dem Rechner 11 verbundenen Monitor 12 wiedergegeben wer- 
den kann. Im Falle des vorliegenden Ausf uhrungsbeispiels ist 
das Datenerf assungssystem 9 mit einer elektrischen Leitung 8, 

20 die in nicht dargestellter Weise beispielsweise ein Schleif- 
ringsystem oder eine drahtlose Ubertragungsstrecke enthalt, 
mit dem Strahlungsdetektor 4 und mit einer elektrischen Lei- 
tung 10 mit dem Rechner 11 verbunden. 

5 Der in der Figur 1 gezeigte Computertomograph kann sowohl zur 
Sequenzabtastung als auch zur Spiralabtastung eingesetzt wer- 
den . 

Bei der Sequenzabtastung erfolgt eine schichtweise Abtastung 
30 des Patienten 3. Dabei wird die Rontgenstrahlenquelle 1 und 
der Strahlungsdetektor 4 beziiglich der Systemachse 6 um den 
Patienten 3 gedreht und das die Rontgenstrahlenquelle 1 und 
den Strahlungsdetektor 4 umfassende Messsystem nimmt eine 
Vielzahl von Projektionen auf, um eine zweidimensionale 
35 Schicht des Patienten 3 abzutasten. Aus den dabei gewonnen 
Messwerten wird ein die abgetastete Schicht darstellendes 
Schnittbild rekonstruiert . Zwischen der Abtastung aufeinan- 
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derf olgender Schichten wird der Patient 3 jeweils entlang der 
Systemachse 6 bewegt. Dieser Vorgang wiederholt sich so lan- 
ge, bis alle interessierenden Schichten erfasst sind. 

Wahrend der Spiralabtastung dreht sich das die Rontgenstrah- 
lenquelle 1 und den Strahlungsdetektor 4 umfassende Messsys- 
tem beziiglich der Systemachse 6 und der Tisch 7 bewegt sich 
kontinuierlich in Richtung des Pfeils b, d.h. das die Ront- 
genstrahlenquelle 1 und den Strahlungsdetektor 4 umfassende 
Messsystem bewegt sich relativ zum Patienten 3 kontinuierlich 
auf einer Spiralbahn c, so lange, bis der interessierende Be- 
reich des Patienten 3 vollstandig erfasst ist. Dabei wird ein 
Volumendatensatz generiert, der im Falle des vorliegenden 
Ausf uhrungsbeispiels nach dem in der Medizintechnik ublichen 
DICOM-Standard kodiert ist. 

Im Falle des vorliegenden Ausf uhrungsbeispiels wird mit dem 
in der Figur 1 dargestellten Computertomographen ein aus meh- 
reren auf einander f olgenden Schnittbildern bestehender Volu- 
mendatensatz 20 des Bauchraums des Patienten 3 angefertigt. 
Der Volumendatensatz 20, der in der Figur 2 schematisch dar- 
gestellt ist, umfasst im Falle des vorliegenden Ausfuhrungs- 
beispiel ca. 250 CT-Schichten (Schnittbilder ) der Matrix 
512x512. In der Figur 2 sind exemplarisch sieben Schnittbil- 
der, die mit den Bezugszeichen 21 bis 27 versehen sind, ange- 
deutet . 

Im Falle des vorliegenden Ausf uhrungsbeispiels sollen abge- 
bildete korperoberf lachennahe Strukturen, die mit dem Volu- 
mendatensatz 20 abgebildet sind, herausgef iltert werden, so 
dass moglichst nur abgebildete innere Organe und die abgebil- 
dete Wirbelsaule des Patienten 3 sichtbar werden. Dazu lauft 
im Falle des vorliegenden Ausf uhrungsbeispiels auf dem Rech- 
ner 11 ein geeignetes Rechnerprogramm, das die nun beschrie- 
bene Schritte ausfuhrt. 
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Zunachst wird in einem ersten Durchgang zur Ermittlung der 
abgebildeten Korperoberf lache jedes Schnittbild 21 bis 27 des 
Volumendatensat zes 20 nach Polarkoordinaten beziiglich einer 
Geraden G, die durch das dreidimensionale Abbild des Bauch- 
5 raums des Patienten 3 verlauft, transf ormiert . Die Gerade G 
ist wenigstens im Wesentlichen rechtwinklig zu den einzelnen 
Schnittbildern 21 bis 27 ausgerichtet . Die Gerade G verlauft 
im Falle des vorliegenden Ausf iihrungsbeispiels im Wesentli- 
chen durch das Zentrum des Volumendatensat zes 20 und ent- 
10 spricht der Z-Achse des den Volumendatensat z def inierenden 
Koordinatensys terns K . 

Jedes Schnittbild 21 bis 27, von denen das Schnittbild 21 in 
der Figur 3 exemplarisch dargestellt ist, ist im Falle des 
15 vorliegenden Ausf iihrungsbeispiels mit kartesischen Koordina- 
ten (x, y) beschrieben. Anschliefiend werden die Bildinforma- 
tionen jedes Schnittbildes 21 bis 27 radial neu angeordnet, 
indem sie beziiglich der Gerade G bzw. beziiglich der jeweili- 
gen Schnittpunkte zwischen der Geraden G und dem entsprechen- 
20 den Schnittbild nach Polarkoordinaten (r, (()) transf ormiert 

werden. Als Beispiel ist der Schnittpunkt S zwischen der Ge- 
raden G und dem Schnittbild 21 in der Figur 3 dargestellt. 
Mit der Transformation nach Polarkoordinaten (r, wird auch 
das Abbild der Korperoberf lache des Patienten 3 transf ormiert 
und als Kontur in jeder transf ormierten axialen Schicht 
(Schnittbild) dargestellt. Eine dem Abbild der Korperoberf la- 
che des Patienten 3 zugeordnete Kontur 40 ist exemplarisch in 
der Figur 4 fur das nach Polarkoordinaten (r, <|>) transfor- 
mierte Schnittbild 21 dargestellt. Das transf ormierte 
30 Schnittbild des Schnittbildes 21 ist mit dem Bezugszeichen 41 
versehen . 



Das Ergebnis der Transformation nach Polarkoordinaten (r, <|>) 
ist ein linear auf getragenes radiales Helligkeit sprof il . In 
35 dieser Rechteckmatrix (abgeleitete Bildmatrix) wird nun eine 
Filterung durchgef iihrt , welche die der Korperoberf lache zuge- 
ordnete Konturen, wie der in der Figur 4 gezeigten Kontur 40, 
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betont. Die Filterantworten ersetzen die Helligkeitswerte in 
der abgeleiteten Bildmatrix. Nun erfolgt die Suche des opti- 
mal en Pfades in dieser Bildmatrix von oben nach unten zum i- 
dentischen Start /Zielpunkt . Das geschieht im Falle des vor- 
liegenden Ausf uhrungsbeispiels mittels dynamischer Optimie- 
rung, wie z.B. in R. Bellman, "Dynamic programming and sto- 
chastic control processes", Information and Control, 1(3), 
Seiten 228-239, September 1958 beschrieben. Der optimierte 
Pfad stellt die radialen Vektoren zu den Korperoberf lachen- 
bildpunkten dar. In einem weiteren Schritt erfolgt eine Ruck- 
transformation der nach Polarkoordinaten transf ormierten Kon- 
turen 40 in die urspriinglichen Koordinaten (x, y, z) des Vo- 
lumendatensat zes 20, so dass das gesamte, durch die einzelnen 
Konturen der Schnittbilder bestimmte Konturensemble und die 
entsprechenden Bildpunkte des urspriinglichen Volumendatensat- 
zes 20 im Zusammenhang der Einzelkonturen iiber alle Schnitt- 
bilder 21 bis 27 uberpriift werden. Dies tragt insbesondere 
zur Unterdruckung von Fehlern (Ausreiftern) und zur Zuverlas- 
sigkeit bei. An vermutlichen Fehlerstellen wird im Falle des 
vorliegenden Ausf uhrungsbeispiels eine Re-Segment ierung in 
den einzelnen Schnittbildern 21 bis 27 mit anschliefiender er- 
neuter Uberprufung des 3D-Kontexts durchgef uhrt . Somit ist 
das Abbild der Korperoberf lache des Patienten 3 im Volumenda- 
tensatz 20 segmentiert. 

Danach erfolgt eine Re-Extraktion rechtwinklig zum Abbild der 
segmentierten Korperoberf lache im Volumendatensat z 20. Wah- 
rend bei der Transformation nach Polarkoordinaten (r, (f>) Hel- 
ligkeitsprof ile rechtwinklig zu alien Punkten eines Kreises 
(idealisierte Oberf lachenkontur ) aus den Originaldaten ermit- 
telt und als Rechteckmatrix aufgetragen wurden, gewinnt man 
bei der Re-Extraktion Profile rechtwinklig zu jedem Bildpunkt 
des Abbildes der segmentierten Korperoberoberf lache (Korper- 
oberf lachenkontur ) . Diese Re-Extraktion wird erneut als 
Rechteckmatrix aufgetragen. Dadurch wird der Volumendatensat z 
20 derart transf ormiert, dass das segmentierte Abbild der 
Korperoberf lache des Patienten 3 in eine Ebene transf ormiert 
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wird. Es entsteht ein in der Figur 5 dargestellter Volumenda- 
tensatz 50, der die Struktur eines Voxelquaders hat. Die in 
die Ebene transf ormierte abgebildete Korperoberf lache ist in 
der Figur 5 mit dem Bezugszeichen 51 versehen und wird nach- 
5 folgend als Mittelebene 51 bezeichnet. Besteht der Volumenda- 
tensatz 20 aus horizontalen Schnitten, wie den Schnittbildern 
21-27, so entspricht jede senkrechte Linie in der Mittelebene 
51 des re-extrahierten Volumendatensat zes 50 ( rechtwinkliger 
Voxelquader) den Bildpunkten des Abbildes der Korperoberf la- 

10 che in je einem der Schnittbilder 21-27. Links von dieser 
Mittelebene 51 im Bereich hoherer y f Koordinaten befinden 
^ sich die CT-Messwerte in der Nahe der Korperoberf lache nach 

* innen. In dem die Form eines Voxelquaders aufweisenden Volu- 
mendatensatz 50 ist der 3D-Zusammenhang fur eine konsistente 

15 Segmentierung gewahrleistet . Er eignet sich daher relativ gut 
fur eine Filterung, urn nichtinteressierende, im Volumendaten- 
satz 50 abgebildete Strukturen herauszuf iltern . Im Falle des 
vorliegenden Ausf iihrungsbeispiels handelt es sich bei den ab- 
gebildeten interessierenden Strukturen urn im Volumendatensat z 

20 50 abgebildete innere Organe und die abgebildete Wirbelsaule 
52. 



Fur das Herausf iltern der im Volumendatensat z 50 nichtinte- 
ressierenden abgebildeten Strukturen werden im Falle des vor- 
5 liegenden Ausf iihrungsbeispiels verschiedene Filteroperat ionen 
fur unterschiedliche nichtinteressierende abgebildete Struk- 
turen ermittelt. Die Filteroperationen berucksicht igen u.a. 
spezifische, den einzelnen nichtinteressierenden Strukturen 
zugeordnete Merkmale und entsprechende Abstandsgewichtungen 

30 zur Korperoberf lache des Patienten 3. Fur tief erliegende Or- 
gane, deren Oberflachen sich der Korperoberf lache annahern, 
werden Mer kmalsf ilterungen angewandt, welche die Wahrschein- 
lichkeit der Zugehorigkeit zur nichtinteressierenden Gewebe- 
schicht verringern. Auch hier wird der Abstand zur Korper- 

35 oberf lache gewichtend mit beriicksichtigt . Zu den Merkmalen 

gehbren auch solche, die mittels dif f erenzierender Operatoren 
ermittelt werden, zur Erkennung von Gewebekonturen wie z.B. 
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Rippenoberf lachen nach innen Oder Organoberf lachen von innen 
nach aufien. Alle Filteroperationen werden anschlieftend in ei- 
ner Wahrscheinlichkeitsmatrix zusammengef asst . Das wird je 
nach Wirkung additiv oder multiplikativ mit passender Skalie- 
5 rung und geeigneten Gewichten durchgef uhrt . 

Wenn, wie es im Falle des vorliegenden Ausf uhrungsbeispiels 

vorgesehen ist, der Volumendatensat z 50 mit kartesischen Ko- 

ordinaten (x ! , y*, z 1 ) beschrieben ist und fur die Bildpunkte 
10 der abgebildeten Korperoberf lache 51 y'=konst. gilt, d.h. 

dass die abgebildete Korperoberf lache 51 in eine Ebene trans- 
^ formiert ist, so sind die x'-z'-Ebenen des Volumendatensat zes 
\ 50 parallel zur abgebildeten Korperoberf lache 51 ausgerich- 

tet. Die Filterung des Volumendatensat zes 50 und somit die 
15 Merkmalsanderung von interessierenden Strukturen zu nichtin- 

teressierenden (herauszuf ilternden) Strukturen erfolgt daher 

in im Wesentlichen y f -Richtung. 

In der zusammengeset zten Wahrscheinlichkeitsmatrix wird wie 
20 bei der Ermittlung der abgebildeten Korperoberf lache im Falle 
des vorliegenden Ausf uhrungsbeispiels die dynamische Optimie- 
rung verwendet, urn einen optimierten Pfad zwischen der abge- 
bildeten Korperoberf lache 51 und den abgebildeten interessie- 
renden Strukturen, also im Falle des vorliegenden Ausftih- 
(4K5 rungsbeispiels die abgebildeten inneren Organe und die abge- 
bildete Wirbelsaule 52, zu finden. Dies erfolgt im Falle des 
vorliegenden Ausf uhrungsbeispiels in dem Volumendatensat z 50 
zugeordneten Schichtbildern mit einem nachf olgenden Durchgang 
zur Gewahrleistung des Zusammenhangs uber das gesamte Volu- 
30 men. Die Optimierung in den einzelnen Schichtbildern ist ei- 
gentlich eindimensional und daher relativ effizient und 
schnell. Die Herstellung des Zusammenhangs erfolgt in der 
Flachendimension; letztendlich wird jedoch 3D-Inf ormation ge- 
liefert. Insgesamt entsteht eine gemeinsame, im Inneren des 
35 abgebildeten Korperinneren des Patienten 3 liegende 3D-Ober- 
flache der interessierenden Strukturen. Diese 3D-Oberf lache 
ist im Volumendatensatz 50, der die Form eines Voxelquaders 
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aufweist, mit dem Bezugszeichen 53 versehen. Das abgebildete 
Gebiet zwischen der 3D-Oberf lache 53 und der abgebildeten 
Korperoberf lache 51 umfasst die nichtinteressierenden Struk- 
turen. Dieses Gebiet wird anschlieftend aus dem Volumendaten- 
satz 50 entfernt; es entsteht ein weiterer, in der Figur 6 
dargestellter Volumendatensatz 60. 

Anschlieftend werden im Falle des vorliegenden Ausf lihrungsbei- 
spiels die Bildpunkte des moglichst nur die abgebildeten in- 
teressierenden Strukturen 52 umfassenden Volumendatensat zes 
60 in das dem ursprunglichen, die Schnittbilder 21 bis 27 um- 
fassenden Volumendatensatz 20 zugeordnete Koordinatensystem K 
zuriick transf ormiert . Es entsteht also aus dem in der Figur 2 
dargestellten Volumendatensatz 20 ein Volumendatensatz, aus 
dem moglichst viele Bildpunkte, die nichtinteressierenden ab- 
gebildeten Strukturen im Volumendatensatz 20 zugeordnet sind, 
herausgef iltert sind und moglichst die Bildpunkte, die den 
interessierenden Strukturen, also den im Volumendatensatz 20 
abgebildeten inneren Organen und der abgebildeten Wirbelsau- 
le, zugeordnet sind, umfasst. In diesem Volumendatensatz kann 
anschliefiend noch ein Volumenrendering der abgebildeten inne- 
ren Organe durchgefuhrt werden. Das Ergebnis ist ein Bild 80, 
das beispielhaft in der Figur 8 wiedergegeben ist. 

Zum Vergleich zeigt die Figur 7 ein Bild 70, das dadurch ent- 
stand, dass beim ursprunglichen Volumendatensatz 20 ein kom- 
plettes Volumenrendering ohne die hier vorgeschlagene Verar- 
beitung durchgefuhrt wurde. 

Der Vergleich Bilder 80 und 70 verdeutlicht den Vorteil des 
erf indungsgemaBen Verf ahrens : 

- Direkte Darstellung "vorderer innerer" Organe 

- Erheblich klarerer Blick auf innere Organe "in der zweiten 
Reihe" 

- Feinere CT-Wert-Auf losung in den Organen 
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- Praxisgerechte Einstellung von verschiedenen 

(Farb) Abschnitten in der Transf erf unktion, d.h. auch bei 
relativ grober Einstellung werden verschieden 
colorierte Objekte klar getrennt und gut dargestellt, was 
ohne die beschriebene Vorverarbeitung in vielen Umgebungen 
die Verwendung des Volumen Renderns unmoglich machen 
wiirde . 

Der Volumendatensatz wird im Falle des vorliegenden Ausfuh- 
rungsbeispiels mit einem Computertomographen hergestellt und 
liegt in Form mehrerer auf einanderf olgender computerto- 
mographischer Schnittbilder vor. Der Volumendatensatz kann 
aber auch mit anderen bildgebenden Geraten, wie insbesondere 
mit einem Magnetresonanzgerat, einem Rontgengerat, einem Ult- 
raschallgerat oder einem PET-Scanner hergestellt werden. Der 
Volumendatensatz muss auch nicht in Form mehrerer aufeinan- 
derf olgender computertomographischer Schnittbilder vorliegen. 

Der Volumendatensatz 20 umfasst im Falle des vorliegenden 
Ausf uhrungsbeispiels einen Teil des abgebildeten Korpers des 
Lebewesens 3. Der Volumendatensatz kann aber auch abgebildete 
korperfremde Objekte, wie z.B. ein Abbild des Tisches 7 des 
Computertomographen, ein Abbild der Kleidung des Patienten 3 
oder ein Abbild von an dem Patienten 3 und in der Figur 1 
nicht dargestellten Instrumenten umfassen. 

Das erf indungsgemafte Verfahren kann auch fur abgebildete 
technische Objekte verwendete werden. Umfasst das technische 
Objekt z.B. eine Ummantelung oder eine Isolierung, so konnen 
diese als nichtinteressierende Strukturen angesehen werden. 

Das Ausf uhrungsbeispiel hat ebenfalls nur exemplarischen Cha- 
rakter . 



200214682 



18 

Patentanspruche 



1. Verfahren zum Herstellen eines Volumendatensat zes, aufwei- 
send folgende Verf ahrensschritte : 

- Segmentieren der abgebildeten Oberflache eines in einem 
ersten Volumendatensat z (20) abgebildeten ersten Objektes 
(3) , 



10 - Transf ormieren des ersten Volumendatensat zes (20) in einen 
zweiten Volumendatensat z (50) derart, dass die segmentierte 
k ^ abgebildete Oberflache in eine Ebene (51) transf ormiert 

v wird, 

15 - Herstellen eines dritten Volumendatensat zes (60), indem der 
zweite Volumendatensat z (50) derart gefiltert wird, dass 
nichtinteressierende, im zweiten Volumendatensat z (50) ab- 
gebildete Strukturen des ersten Objektes (3) aufgrund 
nichtinteressierenden Strukturen zugeordneter Merkmale und 
20 aufgrund der von der Oberflache zu erwartenden Entfernungen 

der nichtinteressierenden Strukturen herausgef iltert werden 
und interessierende, im zweiten Volumendatensat z (50) abge- 
bildete Strukturen (52) des abgebildeten ersten Objektes 
(3) aufgrund interessierenden Strukturen zugeordneter Merk- 
fflK5 male und aufgrund den von der Oberflache zu erwartenden 

Entfernungen der interessierenden Strukturen erhalten blei- 
ben . 



2. Verfahren nach Anspruch 1, bei dem der erste Volumendaten- 
30 satz (20) wenigstens ein abgebildetes zweites Objekt, das au- 

fterhalb des ersten Objekts (3) angeordnet ist, umfasst und 
das abgebildete zweite Objekt aus dem zweiten Volumendaten- 
satz (50) mit den nichtinteressierenden abgebildeten Struktu- 
ren herausgef iltert wird. 

35 

3. Verfahren nach Anspruch 1 oder 2, bei dem bei der Filte- 
rung des zweiten Volumendatensat zes (50) eine den nichtinte- 
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ressierenden Strukturen und/oder den interessierenden Struk- 
turen zugeordnete dichteorientierten, texturorientierte, kan- 
tenbetonte und/oder morphologische Filterung verwendet wird. 

4. Verfahren nach einem der Anspriiche 1 bis 3, bei dem der 
erste Volumendatensat z (20) in Form mehrerer auf einanderf ol- 
gender computertomographischer Schnittbilder (21-27) vorliegt 
oder als Schicht stapel betrachtet wird, die Bilddaten jedes 
Schnittbildes (21-27) mit kartesischen Koordinaten (x, y) be- 
schrieben sind und fur die Segmentierung der abgebildeten 
Korperoberf lache (51) folgende Verf ahrensschritte ausgefuhrt 
werden : 

- Durchfuhren einer Koordinatentransf ormation fur jedes 
Schnittbild (21-27) nach Polar koordinaten (r, <j>) bezuglich 
einer Geraden (G) , die durch das abgebildete erste Objekt 
(3) verlauft und wenigstens im Wesentlichen rechtwinklig zu 
den einzelnen Schnittbildern (21-27) ausgerichtet ist, 

- Ermitteln der Konturen (40), die in jedem transf ormierten 
Schnittbild (41) abgebildet und der abgebildeten Oberflache 
zugeordnet sind, 

- Rucktransf ormieren der Bildpunkte der ermittelten Konturen 
(40) in das dem ersten Volumendatensat z (21-27) zugeordnete 
Koordinatensystem (x, y, z) und 

- Re-Extrahieren von Bildpunkten entlang der Konturen (40) 
fur die Darstellung der in die Ebene transf ormierten Ober- 
flache (51) des abgebildeten ersten Objektes (3) . 

5. Verfahren nach einem der Anspriiche 1 bis 4, bei dem ein 
vierter Volumendatensat z hergestellt wird, indem die Bild- 
punkte des dritte Volumendatensat zes in das dem ersten Volu- 
mendatensatz (21-27) zugeordnete Koordinatensystem (x, y, z) 
rucktransf ormiert werden. 
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6. Verfahren nach Anspruch 5, bei dem ein dem vierten Volu- 
mendatensatz zugeordnetes Bild (60) mittels Volume Rendering 
(VR) dargestellt wird. 
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Zusammenf as sung 

Verfahren zum Herstellen eines Volumendatensatzes 

Die Erfinclung betrifft ein Verfahren zum Herstellen eines Vo 
lumendatensat zes . Zunachst wird die abgebildete Oberflache 
eines in einem ersten Volumendatensat z (20) abgebildeten ers 
ten Objektes (3) segment iert . Danach wird der erste Volumen- 
datensatz (20) in einen zweiten Volumendatensat z (50) derart 
transf ormiert , dass die segmentierte abgebildete Oberflache 
in eine Ebene (51) transf ormiert wird. Aus dem zweiten Volu- 
mendatensat z (50) wird ein dritter Volumendatensat z (60) her 
gestellt, indem der zweite Volumendatensat z (50) derart ge- 
filtert wird, dass nichtinteressierende, im zweiten Volumen- 
datensat z (50) abgebildete Strukturen des abgebildeten Kor- 
pers aufgrund den nichtinteressierenden Strukturen zugeordne 
ten Merkmalen und aufgrund den von der Oberflache zu erwar- 
tenden Entfernungen der nichtinteressierenden Strukturen im 
Korper des ersten Objektes (3) herausgef iltert werden und in 
teressierende, im zweiten Volumendatensat z (50) abgebildete 
Strukturen (52) des abgebildeten Korpers aufgrund den inte- 
ressierenden Strukturen zugeordneten Merkmalen und aufgrund 
den von der Oberflache zu erwartenden Entfernungen der inte- 
ressierenden Strukturen im Korper des ersten Objektes (3) er 
halten bleiben. 



Figur 6 
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